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Introduction 


The  evaluation  of  contour  integrals  with  analytic  integrands  is  often  ac¬ 
complished  by  moving  the  contour  to  an  equivalent  one,  such  that  the  integrand  is 
better  behaved  or  can  be  approximated  more  easily.  In  particular,  movement  of  the 
original  contour  to  one  that  takes  advantage  of  saddle  points  of  the  integrand  or 
paths  of  descent  or  steepest  descent  is  a  very  fruitful  procedure.  It  is  presumed  that 
the  reader  is  familiar  with  this  technique;  see  references  1-3,  for  example.  However, 
one  of  the  difficulties  of  this  procedure  is  determining  the  locations  of  the  steepest 
descent  contours  (reference  1,  p.  263).  For  complicated  integrands,  especially  those 
involving  branches  of  multivalued  functions,  exact  determination  of  steepest 
descent  contours  is  virtually  impossible  analytically,  and  recourse  to  some  type  of 
computer  aid  is  recommended.  The  procedure  given  here  does  not  require  solution 
of  nonlinear  equations,  but  does  give  a  very  good  indication  of  steepest  descents 
with  a  minimum  of  analytical  and  programming  effort. 
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Explanation  of  Technique 

Suppose  we  wish  to  evaluate  the  contour  integral 

fdz  g ( z )  exp(Xw(z) )  ,  (]) 

where  g(z)  and  w(z)  are  analytic  functions  of  z,  except  for  isolated  singularities  such 
as  poles,  essential  singularities,  and  branch  points;  and  C  is  a  contour  (finite  or 
infinite)  in  the  complex  z-plane.  The  saddle  points  of  the  exponential  in  (1)  occur 
where  (reference  1,  p.  258) 


w'(zs)  =  0  .  (2) 

For  A  real,  the  standard  method  of  determining  the  paths  of  steepest  descent  out  of  a 
saddle  point  is  to  keep  the  imaginary  part  of  w(z)  constant  and  equal  to  its  value  at 
zs  (reference  1,  p.  255).  This  generally  leads  to  difficult  transcendental  equations 
that  must  be  analytically  investigated  approximately  or  solved  numerically. 

An  alternative  procedure  for  finding  the  steepest  descent  directions  at  any  point 
in  the  z-plane  is  as  follows:  Let 

z  =  x  +  iy  , 

,  .  *  (3) 

w(x  +  iy)  =  u  +  iv 

Then  the  magnitude  of  the  exponential  in  (1)  is  exp(Au),  and  its  direction  of  steepest 
descent  at  x,  y  is  proportional  to  the  negative  of  the  gradient  (reference  1 ,  p.  254): 

V  exp (Xu)  =  X  exp (Xu)  Vu  ,  (4) 


where 


Vu 


3U  -► 
3x  ax 


+ 


dU  -* 

9y  ay  * 


(5) 


and!?x  and  a^  are  unit  vectors  in  the  positive  x-  and  y-directions,  respectively.  The 
explicit  evaluation  of  (5)  requires  that  one  analytically  evaluate  u  =  Re{w(x  +  iy)} 
and  then  analytically  derive  3u/3x  and  3u/3y.  This  can  be  tedious  and  is  liable  to 
human  error. 


An  alternative  simpler  procedure  is  possible:  By  the  Cauchy-Riemann  conditions 
applied  to  a  function  analytic  at  z,  the  derivative 


w'(z) 


3u  .  3v 
3x  +  1  3x 


3v  .  _3u 
3y  1  3y 


(6) 


That  is,  we  can  express  the  desired  partial  derivatives  as 


fy  =  Re{w'(z)},  =  -  Im(w'(z)} 


(7) 
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Then  the  negative  of  the  gradient  in  (5)  can  be  written  as 

-Vu  =  -Reiw'fz)}  a  +  Im{w'(z)}  a  .  (8) 

x  y 

So  if  we  evaluate  w  '(z),  the  direction  of  steepest  descent  at  any  point  z  has 
components  in  the  x,  y  directions  proportional  to 

-Re{w'(z)},  Iid{w'(z)  }  (9) 

for  w  '(z)  #  0.  The  only  analytical  calculation  necessary  is  that  of  derivative  w  '(z),  a 
task  generally  easily  accomplished,  and  indeed  necessary  for  evaluation  of  saddle 
point  locations  anyway.  A  computer  program  can  then  be  written  to  numerically 
evaluate  w  '(z)  in  terms  of  its  components  (9)  at  all  points  of  interest  in  the  z-plane. 
A  program  for  this  procedure  is  given  in  appendix  A,  along  with  the  specific 
examples  displayed  later  in  this  report. 

Presentation  of  this  steepest  descent  information  for  human  interpretation  is 
accomplished  here  by  drawing  a  short  standard-size  line  through  each  point 
z  =  x  +  iy,  centered  on  the  point  and  with  an  arrowhead  pointing  in  the  direction 
of  steepest  descent.  The  magnitude  of  the  rate  of  steepest  descent  is  discarded;  only 
the  direction  is  preserved.  How  effective  this  procedure  is  will  be  demonstrated  by 
the  following  examples. 
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Examples 

Airy  Function 

The  Airy  function,  Ai,  is  proportional  to  (1)  when  g(z)  =  1, 

w(z)  =  z  -  ,  (10) 

and  C  is  an  infinite  contour  starting  anywhere  in  the  angular  range  -5n/6  <  arg  (z) 
<-n/2  and  ending  in  n/2  <  arg(z)  <5n/6;  see  reference  1 ,  pp.  52  and  266.  Then 

w'(z)  =  1  -  Z2  .  (11) 

A  computer  need  only  evaluate  the  complex  product  z*z  and  subtract  it  from  1  in 
order  for  (11)  to  be  used  in  (9).  An  example  of  this  procedure  is  given  in  figure  1. 
The  arrows  clearly  indicate  the  steepest  descent  paths  from  any  point  in  the  z-plane. 
The  two  solutions  of  (1 1)  equal  to  zero,  namely,  saddle  points  zs  =  ±1,  have  arrows 
pointing  both  inward  and  outward  at  these  points,  reflecting  the  very  nature  of  a 
saddle  point.  Movement  of  the  original  contour  C  to  the  steepest  descent  contours 
(solid  lines)  out  of  the  saddle  point  at  zs  =  -1  is  easily  accomplished;  no 
singularities  of  the  integrand  of  (1)  are  crossed  in  the  movement  process. 

The  Airy  function  for  complex  argument  (reference  1,  section  7.3)  has,  more 
generally, 

w(z)  =  zexp(iO)-  jz3  .  (12) 

Figure  1  corresponded  to  0  =  0.  For  0  #  0,  the  determination  of  steepest  descent 
paths  is  analytically  difficult  (ibid.).  However,  since 

w'(z)  =  exp(i9)  -  z“  ,  (13) 

computer  evaluation  of  (13)  and  (9)  is  trivial.  The  steepest  descent  directions  for 
0  =  3n/4,  for  example,  are  depicted  in  figure  2.  The  solid  lines  in  the  neighborhood 
of  the  saddle  points  were  hand-drawn  upon  observation  of  the  descent  arrows.  The 
interpretation  of  figure  2  is  much  easier  than  its  counterpart  in  reference  1,  figure 
7.3.3.  Also  the  determination  of  an  equivalent  contour  to  C  is  easily  achieved  by 
reference  to  figure  2.  First  let  the  new  contour  come  from  00  exp(-i2n/3)  along  a 
steepest  ascent  to  the  saddle  point  at  zs  =  exp(-i5n/8).  Then  let  it  continue  in  a 
northeasterly  direction  along  the  steepest  descent  direction  out  of  this  saddle  point 
to  the  y-axis.  Next  proceed  due  north  to  a  point  near  z  =  i  and  then  follow  a 
steepest  descent  path  out  to  °°exp(i2«/3).  The  vertical  portion  of  this  new  contour  is 
not  a  path  of  steepest  descent,  but  it  is  obviously  a  path  of  descent  because  the 
projections  of  the  arrows  on  this  vertical  section  all  point  in  the  upward  direction  of 
travel.  This  discussion  also  points  out  that  the  other  saddle  point  at  zs  =  exp  (i3n/8) 
does  not  enter  into  the  asymptotic  development  of  Ai,  at  least  for  this  value  of  0. 

An  alternative  equivalent  contour  to  C  is  the  pair  of  steepest  descent  contours 
passing  through  the  saddle  points  and  connecting  °°exp(-i2n/3)  to  +°°  and  +°°  to 
°°exp  (i2n/3),  respectively;  see  figure  2.  Observe  that  the  movement  of  C  to  this  new 
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pair  of  paths  is  not  an  approximation;  it  is  an  alternative  exaet  representation  of  the 
original  integral.  If  one  now  approximates  these  two  path  integrals  by  their  con¬ 
tributions  near  their  peaks  at  the  saddle  points,  the  saddle  point  at  zs  =  exp(i3n/8) 
will  yield  exponentially  small  contributions  relative  to  that  at  zs  =  exp(-i5n/8). 
This  can  be  seen  from  figure  2  by  drawing  a  straight  line  between  the  two  saddle 
points;  all  projections  of  arrows  along  this  line  point  at  the  upper-right  saddle  point, 
meaning  that  the  value  of  |exp(Aw(z))|  is  smaller  there. 


Figure  1.  Steepest  Descent  Directions  for  Airy  Function,  6  =  0 
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Figure  2.  Steepest  Descent  Directions  for  Airy  Function,  9  =  3n/4 


When  we  continue  on  to  the  case  8  =  n  in  (12),  the  character  of  steepest  descents 
is  as  depicted  in  figure  3.  Now  the  straight  line  connecting  the  two  saddle  points  has 
all  arrows  perpendicular  to  it;  thus  this  vertical  line  is  a  contour  of  constant 
|exp(A\v(z))|.  This  means  that  both  saddle  points  contribute  equally  to  the  value  of 
integral  (1).  Again  the  pair  of  steepest  descent  contours  through  the  saddle  points 
(mentioned  in  the  above  paragraph)  represent  exactly  the  original  integral;  one 
could  evaluate  the  original  integral  exactly  by  adding  the  total  contributions  of  both 
of  these  paths,  or  an  approximation  can  be  achieved  by  computing  the  integrand 
near  its  peaks  at  the  saddle  points. 
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Figure  3.  Steepest  Descent  Directions  for  Airy  Function,  6  =  n 
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Hankel  Function 

The  Hankel  function  takes  the  form  (reference  1, 7.2.23) 

w(:)  =  i(cos(z)  +  .  (z  -  T)J  ,  (14) 

where  ft  is  a  constant.  There  follows  immediately 

w'(z)  =  i [ -  sin (z)  +  3]  .  (15) 

Computer  evaluation  of  (15)  requires  only  a  trigonometric  sin  of  a  complex  number, 
followed  by  subtraction  and  multiplication.  Sample  descent  direction  plots  for  ft  = 
0.5,  1,  and  1.5  are  given  in  figures  4-6,  respectively.  (It  is  informative  to  compare 
these  figures  with  figures  7.2. 1-7. 2. 4  in  reference  1.)  Movement  to  equivalent 
contours  is  obvious  from  figures  4-6.  Since  (15)  has  period  2n  in  x,  only  a  2n  strip 
has  been  plotted  in  figures  4-6.  The  character  of  the  steepest  descents  in  figure  5  for 
P  =  1  is  different,  in  that  the  saddle  points  have  coalesced;  however,  there  is  no 
difficulty  ascertaining  from  the  plots  what  new  contour  to  adopt. 

Klein-Gordon  Equation 

This  example  is  complicated  by  the  presence  of  branch  lines  in  the  exponent 
function  w(z).  Specifically  we  have  (reference  1 , 7.5.9) 

w(z)  =  i[e(z2  -  nj/2  -  z]  ,  (16) 

where  the  branch  of  the  square  root  is  taken  as  positive  real  for  z  =  x  >  1 ,  and  with 
branch  lines  extending  vertically  downward  from  the  branch  points  at  z  =  ±1.  The 
derivative  of  (16)  is 


where  the  same  square  root  branch  must  be  taken  as  in  (16). 

If  one  has  available  a  computer  program  that  evaluates  the  principal  square  root 
of  a  complex  number,  denoted  here  by  zl/;,  it  can  be  used  to  evaluate  (17)  in  the 
following  manner.  Observe  first  that  the  branch  line  of  principal  square  root  z'- 
occurs  where  z  =  -p  for  p  ^  0;  i.e.,  p  can  take  on  all  nonnegative  real  values.  So 
consider  the  representation 

(z2  -  l)^2  =  i (-iz  +  i)1/2  C-iz  -  i)1/2  .  (18) 

For  z  real,  positive,  and  large,  the  right-hand  side  of  (18)  approaches 
i[exp(-in/4)z';]  [exp(-in/4)z  ;]  =  z,  as  desired.  Furthermore,  the  two  branch  lines 
of  (18)  occur  where  the  arguments  of  the  two  principal  square  roots  have  values 

-iz  +  i  =  -p  for  p  >_  0  ; 

i.e.,  z  =  +.1  -  ip  for  p  ^  0  .  (19) 
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Figure  4.  Steepest  Descent  Directions  for  Hankel  Function,  /? 
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Figure  6.  Steepest  Descent  Directions  for  Hankel  Function,  fi  =  1.5 
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These  are  vertieally  downward  from  z  =  ±1,  as  desired.  Thus  (17)  can  be  easily 
evaluated  by  taking  the  two  complex  principal  square  roots  indicated  in  (18)  and 
performing  multiplication,  division,  and  addition  of  complex  numbers.  The  specific 
coding  is  illustrated  in  appendix  A. 

The  pole  of  g(z)  at  z  =  v0  for  this  example  (reference  1,  7.5.8)  has  no  effect  on 
the  steepest  descent  contours  of  w(z).  The  steepest  descent  directions  for  6  =  0.8 
are  depicted  in  figure  7.  There  are  saddle  points  at  zs  =  ±5/3,  and  the  steepest 
descent  contours  go  vertically  downward  eventually.  The  steepest  descent  directions 
near  the  branch  lines  emanate  from  the  branch  lines  themselves,  but  these  branch 
lines  have  no  effect  on  the  steepest  descent  contours  through  the  saddle  points. 
However,  if  the  branch  line  emanating  from  the  branch  point  at  z  =  1  had  been 
taken  at  angle  -n/6,  for  example,  it  would  have  interfered  with  the  steepest  descent 
contours  in  the  4th  quadrant  of  the  z-plane.  Such  a  choice  of  branch  for  the  square 
root  in  (16)  and  (17)  is  undesirable  and  should  be  avoided,  as  was  done  in  figure  7. 

Function  with  Essential  Singularity 

This  integrand  is  characterized  by  (reference  4) 


exp(w(z))  =  z  exp 


(i4t)  ■ 


which  function  has  a  zero  at  z  =  0  and  an  essential  singularity  at  z  =  l.Then 

w(z)  =  In (z) -  -■  ■ K  -j-  , 


,  1  K 

w  (z)  +  - - -y  .  (21) 

(z  -  1) 

The  steepest  descent  contours  for  (21)  are  depicted  in  figure  8  for  K  =  3.  The 
essential  singularity  generates  a  “dipole effect”  about  z  =  l,i.e.,0atz  =  l  +  ,and 
°°  at  z  =  1-,  for  K  >  0.  The  zero  of  (20)  at  z  =  0  manifests  itself  as  a  point  toward 
which  all  the  arrows  point,  since  zero  is  the  smallest  magnitude  that  any  complex 
function  can  take  on.  The  saddle  points  (roots  of  (21))  occur  at 

zg  =  exp(+.i9),  where  9  =  arc  cos^l  -  -yj  ,  for  0  K  4  .  (22) 

Utilizing  the  information  in  figure  8,  we  find  it  relatively  easy  to  decide  what  the 
effect  of  moving  an  original  contour  around  in  the  z-plane  will  do  to  integral  (1). 
Movement  across  the  essential  singularity  at  z  =  1  will  necessitate  consideration  of 
the  residue  at  this  point.  The  zero  and  saddle  points  of  (20)  are  not  points  of 
singularity. 

Cubic  Function 

This  example  comes  from  reference  3,  pp.  296-302;  it  is  characterized  by 
g(z)  =  Vi  and  (ibid.,  upper  line  of  6.6.25) 


w(z)  =  -e^a  ^4z^  -  i5z  +  iz3)  . 
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Figure  8.  Sleepesl  Descent  Directions  for  exp  (w(z))  =  z  exp  (— j  j,  K  =  3 

Also  now  the  integral  is  a  finite  one,  from  z  =  -ltoz  =  +  1.  This  example  exhibits 
a  Stokes  phenomenon  at  certain  values  of  a,  where  we  have  represented  \  =  | A|e'“. 
We  find 

wf(z)  =  -e*a  ^8z  -  i5  +  iSz*-)  ,  (24) 

which  has  zeros  (saddle  points)  at  zs  •=  i  and  i5/3. 

A  plot  of  steepest  descent  directions  for  a  =  0  (positive  real  A)  is  given  in  figure  9. 
It  indicates  that  the  steepest  descent  contours  out  of  the  limits  at  z  =  - 1 ,  +  1  tend  to 
<»exp(-i5n/6)  and  »exp(-in/6),  respectively.  But  these  two  valleys  at  00  can  be 
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Figure  9.  Steepest  Descent  Directions  for  Cubic,  a  =  0 


joined  by  the  saddle  point  contribution  through  the  point  zs  =  i.  Thus  the  integral 
over  (-1,  1)  is  exactly  equal  to  the  sum  of  these  three  steepest  descent  contours.  The 
saddle  point  at  zs  =  i5/3  need  not  be  considered.  The  dominant  contribution  is 
obviously  that  at  the  saddle  point  zs  =  i,  as  may  be  seen  by  the  arrow  directions. 

For  a  =  5n/12,  figure  10  indicates  a  similar  behavior.  Since  the  steepest  descent 
contours  out  of  -1,  +1  tends  to  -175°  and  -55°, respectively,  the  saddle  point  at 
zs  =  i  must  again  be  used  to  join  them.  The  dominant  contribution  is  seen  to  be  due 
toz  =  -1,  by  making  use  of  the  arrows  of  descent  in  this  figure. 

For  a  =  3n/4,  however,  we  see  from  figure  1 1  that  both  of  the  steepest  descent 
contours  out  of  -1,  +1  tends  to  °°exp(-i5n/12).  Now  there  is  no  need  to  employ 
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either  of  the  saddle  points.  The  integral  over  (-1,1)  is  given  exactly  by  the  sum  of  the 
two  steepest  descent  contributions.  The  dominant  contribution  is  again  due  to 
z  =  -1,  since  we  have  to  descend  from  z  =  -1  to  near  z  =  0  to  reach  magnitude 
values  comparable  to  those  at  z  =  +1. 


Figure  10.  Steepest  Descent  Directions  for  Cubic,  a  =  5n/12 


16 


*  *  *  *  t 
\  s  s  , 

•j-  **  y 

4-  ^  -  —  — 
4.  v  «*  ■*- 

{  V  \  V  \ 

i  \  \  \  \ 

4-v-v-v-v 


ii 


(«k  w-  —  ^.-  _ — r — / 

■ - \  *  /  t  f 

^  j  *  /  /  / 

\  s  y  /  f  f  t 

i  —  -*  s  J  t  t  t  t 

4—^  -/*-/  i-f --*-+-+■ 

i^v  t  t  )  t  t  \  \ 

•fSADDLE  \  i  \  \  \  \ 

4  W\  X  4  N  V  V  X 

\vvvv\VV\V 

vL  SADDLE  ! 

T  POINT  1 

V  i  ✓  - 

>V  *  '  '  T - 

4  \  (  /  /  *  ✓ - 

4  V\  I  iff  BTrcnper' 


t 

i 

i 

f 

n  f  f  s  '  \ 

i 

i 

t 

/ 

\  f  f  f  f  4 

* 

l 

i 

/ 

t  f  /  f  f  4 

\ 

l 

i 

f 

\  t  f  f  f  4 
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Gaussian  Exponent 


The  characteristic  function  of  a  particular  type  of  impulsive  noise  is  given  in 
reference  5,  equation  (5),  in  the  form 

°o  m  /  ,  mB“  A 

f(0  =  cxp(-A)  £  ^  exp  I-  J  T?  I  (25) 

m=0  \  / 


for  a  purely  Poisson  process  (no  Gaussian  background).  This  summation  can  be 
evaluated  in  the  closed  form 
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f(0 


(26) 


The  corresponding  probability  density  function  of  this  impulsive  noise  is  given  by 
the  Fourier  transform 


p(v)  =  37  /  d"  CXP 


+  A 


i  1  0 

exP  "  T  — 


\ 

) 


By  expanding  f(4)  in  (26)  in  a  power  series,  w£note  that  the  mean  square  value  of 
the  Poisson  process  is  readily  found  to  be  AB2/2.  Since  we  will  be  interested  in 
values  for  the  dimensionless  parameter  A  of  the  order  of  1  (e.g.,  A  =  0.35  in 
reference_  5,  figure  3),  we  will  normalize  our  random  variable  according  to 
t  =  v/(B12,/2)l/-  =  v/o0.  Also,  as,  4  “*■  ±°°,  f(4)  in  (26)  tends  to  nonzero  value 
exp(-A).  Adding  and  subtracting  this  quantity  and  letting  4  =  z/o0  enables  (27)  to 
be  expressed  as 


p(v)  -  exp(-A)  6(v) 


+  2  Vo~^  /  dZ  exP(_itz)[exp(A  exp(-z2/2) )  -  l]  . 

o  (28) 

Although  there  is  no  obvious  parameter  k  in  this  form,  it  is  shown  in  appendix  B 
that  we  can  still  use  steepest  descent  procedures  on  the  logarithm  of  the  integrand  of 
(28);  i.e.,  here  we  have 

w(z)  =  -itz  +  In [exp(a (z) )  -  1]  ,  (29) 

where  we  have  defined 


a(z)  =  A  exp(-z2/2)  .  (30) 

Observe  that  as  z  -*•  00  with  arg(z)  in  the  two  sectors  within  n/4  of  the  positive-real 
or  negative-real  axes,  a(z)  becomes  very  small  and 

w(z)  ~  -itz  +  ln[a(z)]  =  -itz  +  In  A  -  z2/2  ;  (31) 

this  means  that  Re  w(z)  -*•  in  these  two  sectors  of  the  z-plane,  and,  therefore,  the 
integrand  of  (28)  tends  to  zero  in  these  sectors  as  z  -*  °°. 

The  integrand  of  the  integral  of  interest  in  (28)  is  zero  in  the  finite  z-plane  only 
when 

h(zQ)  =  exp^Ae  -  1  =  0  , 

2  n 

Ae  Z°  =  i2wn  for  n  jl  0  , 
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2 

z 

o 


2 


In  li 


+  i2-rrm 


i—  sgn(n)  +  i2Trm 


z 

o 


+  iir(4m  +  sgn(n)) 


«• 


1/2 


(32) 


where  A  >  0,  the  square  root  is  the  principal  branch,  and  n  and  m  are  arbitrary 
integers  (negative,  zero,  or  positive),  except  that  n  ^  0.  These  zeros  of  h(z)  are 
important  because  they  will  be  locations  toward  which  all  the  steepest  descent 
arrows  must  point  in  their  neighborhoods,  since  zero  is  the  smallest  magnitude  that 
a  complex  function  can  take  on.  These  zero  locations,  of  which  there  are  an  infinite 
number,  depend  only  on  A,  and  not  on  normalized  variable  t  =  v/oQ  in  expression 
(28)  for  the  probability  density  p(v).  The  locations  of  the  zeros  of  h(z)  in  the  third 
quadrant  closest  to  the  origin  are  depicted  in  figure  12  for  A  =  0.35.  There  is 
symmetry  in  the  other  quadrants  since  (from  the  first  line  of  (32)) 


h(-z)  =  h(z),  h(z*)  =  h*  (z)  .  (33) 

For  purposes  of  evaluating  steepest  descent  directions  and  saddle  point  locations, 
we  note  that  a(z)  in  (30)  has  the  property 

a/  (z)  =  -za(z)  ,  (34) 


and  so  (29)  yields 


w’  (z)  =  -it 


_ z  a(z) 

1  -  exp(-a(z)) 


(35) 


Thus  the  saddle  points,  zs,  of  which  there  are  an  infinite  number,  are  solutions  of 


zs  a(zs) 

1  -  exp(-a(zs))  ’  (36) 

and  obviously  depend  on  both  A  and  t.  The  most  important  saddle  point  is  at 
zs  =  -i ft,  P  positive  real,  where  (using  (30)), 

_ BA  exp(B"/2)  _  _ 

1  -  exp(-A  exp(6“V2))  (37) 

The  steepest  descent  directions  for  A  =  0.35,  t  =  2  are  depicted  for  (35)  in  figure 
13.  The  saddle  point  at  zs  =  -ip  =  -il  .341  satisfies  (37);  there  are  five  other  saddle 
points  indicated  by  X  in  the  figure,  in  addition  to  the  zeros  carried  over  from  figure 
12.  The  steepest  descent  contour  out  of  the  saddle  point  zs  =  -il  .341  is  drawn  as  a 
solid  line;  it  is  asymptotic  to  y  =  -t  as  z  -*  <»  with  arg(z)  within  n/4  of  the  positive- 
and  negative-real  axes.  Movement  of  the  original  contour  from  the  real  axis  in  (28) 
to  the  steepest  descent  contour  is  easily  justified. 
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Figure  12.  Zero  Locations  fur  A  =  0.35 
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Figure  13.  Steepest  Descent  Contour  for  A  =  0.35,  t 
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To  deduce  the  asymptotic  nature  of  the  steepest  descent  contours,  recall  (31);  thus 
the  imaginary  part  of  w  is 

v~-tx  -  xy  =  -x(y  +  t)  as  z  (38) 

in  the  sectors  under  consideration.  But,  since  from  (29), 

w(-i.  )  =  -,  t  +  In [exp (A  exp (  '“/-))  -  0  (39) 

is  real,  we  require  (38)  to  approach  0  as  x  -*  ±°°.  This  requires  that  y  -*•  -t,  as 
claimed. 

We  observe  from  figure  13  that  no  use  is  made  of  the  saddle  point  at  (-1.71, 
-2.38),  nor  of  the  infinite  number  of  other  saddle  points.  We  also  observe  that  the 
descent  directions  in  the  neighborhood  of  the  closely  spaced  zeros  and  saddle  points 
near  the  bottom  of  the  figure  is  very  detailed  and  complicated;  however,  none  of 
that  information  is  needed. 

When  t  is  increased  to  5,  figure  14  applies.  Now  the  steepest  descent  contour  out 
of  the  saddle  point  at  zs  =  -i/3  =  -il  .945  heads  into  the  zero  at  (-0.632,  -2.485). 
How  to  connect  from  this  latter  point  to  z  =  is  not  clear.  Instead,  we  consider  a 
horizontal  descent  contour  out  of  zs  =  -il  .945  until  we  get  in  the  neighborhood  of 
(-2,  -2),  and  then  we  resume  a  steepest  descent  contour  heading  toward 
y  =  -t  =  -5.  The  major  contribution  to  this  descent  contour  is  given  by  the 
neighborhood  of  the  saddle  point  at  x  =  0;  the  descent  and  steepest  descent  con¬ 
tours  are  tangent  at  this  saddle  point. 

The  aid  afforded  by  the  steepest  descent  directions  depicted  in  figures  13  and  14  is 
extremely  worthwhile,  since  the  exponential  in  (29)  and  (35)  makes  an  analytical 
approach  very  difficult.  The  ability  to  discard  or  avoid  certain  regions  of  the  z-plane 
in  determining  an  appropriate  descent  contour  is  rather  obvious  from  the  figures 
when  coupled  with  basic  information  about  the  integrand,  like  the  asymptotic 
behavior  of  the  steepest  descent  contour. 

The  above  results  yield  exact  values  for  the  original  integral,  since  we  have  simply 
determined  contours  equivalent  to  those  originally  specified.  Now  we  will  derive  an 
asymptotic  expansion  for  the  probability  density  in  (28)  for  large  t  =  v/o0,  i.e.,  at 
values  v  much  larger  than  the  standard  deviation  of  the  Poisson  process. 

To  do  this,  we  need 


where  a  =  a(z)  is  defined  in  (30).  Then  we  find 

w"(Zs)  =  t2  ♦  if-  (l  -  z2s-  z2s  a(zs))  (41) 

s 

and,  in  particular, 

w"(-i6)  =  t2  -  +  B2  +  B2  A  exp(82/2))  ,  (42) 
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for  the  saddle  point  at  zs  =  -i/3.  Also  w(-i/3)  is  given  by  (39).  Then  we  have  the 
approximation  (reference  1 ,  chapter  7) 


(43) 

One  drawback  with  this  solution  is  that  fi  depends  on  t  through  the  solution  of 
transcendental  equation  (37).  For  large  t,  we  have,  to  first  order, 

8  as  [2  ln(t/A)]1/2  5  L1/2  .  (44) 

For  development  of  additional  terms  and  the  general  philosophy  of  solution  of 
these  types  of  problems,  see  reference  6,  pp.  11-16  and  83-84.  We  find,  more 
generally, 

e2  "  L  -  r^T  In  (L)  .  (45) 

Even  so,  substitution  into  (43)  yields  a  very  complicated  expression  for  the 
probability  density  function  unless  t  is  excessively  large.  It  can  be  seen  that  (43) 
decays  slightly  faster  than  an  exponential  for  large  t. 
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Summary 

A  technique  for  simple  determination  of  the  steepest  descent  direction  at  any 
point  in  the  complex  plane  has  been  presented  and  illustrated  with  numerous 
examples.  The  movement  of  the  original  contour  to  an  equivalent  descent  or 
steepest  descent  contour  is  an  exact  representation  and  can  be  deduced  fairly  easily 
from  the  descent  information.  At  this  point,  two  alternatives  are  available,  either 
exact  numerical  evaluation  of  the  integral  or  an  approximation  such  as  Laplace’s 
method.  Very  difficult  integrands  can  be  handled  very  effectively  via  this  approach. 
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Appendix  A 

Computer  Programs 

The  major  part  of  the  calculations  required  for  the  six  examples  in  the  main  text  is 
that  of  w  '(z).  This  is  accomplished  in  the  Subroutine  Wderivative  in  the  main 
program  written  in  BASIC  and  listed  at  the  end  of  this  appendix.  But,  first,  the  six 
subroutines  for  the  examples  given  are  listed.  They  illustrate  how  little  program¬ 
ming  is  actually  needed  to  compute  w  '(z),  provided  that  one  has  already  written 
subroutines  for  the  standard  complex  operations  and  functions  like  multiply, 
divide,  exp(z),  log(z),  square  root,  arg(z),  sin(z),  cos(z),  etc.  These  latter  functions 
are  listed  for  completeness  as  subroutines  at  the  end  of  the  enclosed  program. 

41®  SUE  Wder  i  vat  i  oe<X,  Y,  Rewl ,  Imwl )  !  Airy  function 

4 1 1  COM  Ct ,  St-  !  c os(  t  bet  a> ,  ■=  i  n  <  t  be  t 

412  CALL  Mu I <X,  Y, X, Y, A, B> 

420  R  e  w  1  =  C  t  -  A 

430  Imwl=St-B 
440  SUBEND 


410  SUB  Wder  i  oat  i  oe  <  X ,  Y ,  Rewl ,  I  mu  1  >  1  Hank  *  1  f  ur.c  t  i  on 

411  Beta*. 5 

412  CALL  Si  ri<K,  Y,  A,  B> 

420  Rewl=B 

430  I  m  w  1  =  B  e  t  a-  A 

440  SUBEHD 


410  SUE  Wder i oat i oe<X, Y, Rewl , Imwl >  !  K  1  e  i n  -Gordon 

411  Theta*. 8 

412  CALL  Sqr <Y, -X+l , ft, B> 

413  CALL  Sqr  < Y , -X- 1 , C , D  > 

414  CALL  Mul <A,B,C,D,E,F> 

415  CALL  Dio<Theta*X,Theta*Y,E,F,i3,H> 

420  Rewl=G 

430  Imwl=H-l 
440  SUBEND 


410  SUB  Hder  i  oat  i  oe  (.  X,  Y ,  Rewl ,  I  mw  1  >  1  Essential  singularity 

411  K  =  3 

412  CALL  Mul < X- 1 , Y, X- 1 , Y, A, B > 

413  CALL  D i o < K , 0 , A , B , C , D ) 

414  CALL  D i o ( 1 , 0 , X , Y , E , F ) 

420  Rewl*C+E 

430  Imwl*D+F 
440  SUBEND 
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SUB  Wde-r  i  vat  i  we  <  X ,  Y ,  Rew  1 ,  I  mw  1  !  cubic  fur ct ion 

COM  Ca,Sa  !  costal  pha),  s  i  ntal  pKa> 

CAUL  Mu  1 <X, Y,X, Y,A,B> 

CALL  Mu  1 < Ca, Sa, 8*X-3*B, 8*Y-5  +  3*A,  C,  D> 

R  e  w  1  =  -  C 
I hh 1 =-D 
SUBEND 


SUB  Wd«r  i vat i ve<X, Y, Rewl , Irnwl > 
A=.  35 
T  =  2 

CALL  Mul <X,Y,X,Y,T1,T2> 

CALL  Exp<-.5*T1,-.5*T2, T3, T4> 

Ar=A*T3 

Ai =A*T4 

CALL  Exp< -Ar , -Ai , Er , E i > 

CALL  Mul<X,Y,Ar,fli ,T1,T2) 

CALL  D i v < T 1 , T 2 , 1-Er,-Ei , T3,T4> 

Rew 1 =-T3 

Imwl=-T-T4 

SUBEND 


!  Gaussian  e*pon€ri t 
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1  ! 

Steepest  Descent  vie  w" ;  use  SUB  Mderiuative  i 

18 

X 1  =  -  3 

!  LEFT  ABSCISSA 

20 

X  2  =  3 

!  RIGHT  ABSCISSA 

30 

Y 1  =  -  3 

'  BOTTOM  ORDINATE 

40 

Y  2*3 

!  TOP  ORDINATE 

50 

B  x  = .  2 

!  X  INCREMENT 

60 

By*.  2 

!  Y  INCREMENT 

70 

PLOTTER  IS  "GRAPHICS" 

88 

GRAPHICS 

90 

SCALE  X 1 , X2 , Y 1 , Y2 

100 

LINE  TYPE  3 

110 

GRIB  1,1 

!  GRID  LINE  SPACING 

120 

LINE  TYPE  1 

. 

130 

F  * . 2*SQR  < Bx*Bx+Dy*Dy ) 

!  ARROW 

140 

G 1  =  1 -COS  <  P I  1 2  >  * .  8 

!  INFOR- 

150 

G2  =  S  I N  ( P I  •'  1 2  >  * .  8 

!  MAT  I  ON 

160 

FOR  X=X1  TO  X2  STEP  Dx 

170 

FOR  Y=Y 1  TO  Y2  STEP  By 

180 

CALL  Wder i vat i oe t X, Y, Rewl , Imwl >  !  w / < z > 

190 

CALL  Di rect i on<Rewl , Imwl , Cos, Si n>  !  direction 

200 

IF  ABS< Cos > +ABS< S i n > >0 

THEN  230 

210 

OUTPUT  0 ;  11 SABBLE  POINT 

AT  “ ; X ; Y 

220 

GOTO  350 

230 

T 1 *F*Cos 

240 

T2*F*Sin 

250 

Xa=X+Tl*Ql 

260 

Xb=T2*Q2 

270 

Y a*Y+T2*Ql 

280 

Yb=Tl*Q2 

290 

MOVE  X-T1.Y-T2 

300 

DRAW  X+Tl , Y+T2 

310 

MOVE  Xa+Xb, Ya-Yb 

320 

DRAW  X+Tl , Y+T2 

330 

DRAW  Xa-Xb,  Ya*-Yb 

340 

PENUP 

350 

NEXT  Y 

360 

NEXT  X 

370 

PAUSE 

380 

DUMP  GRAPHICS 

390 

END 

400 

! 

410 

SUB  W  d  e  r  i  \>  at  i  m  e  <  X ,  Y ,  R  e  w  1 ,  I  m  w  1  )  !  A  i  r  y  f  u  n  c 

420 

R  e  w  1  =  1  -  X  *  X + Y  *  Y 

1  W<.  Z)*2-2  ' 

430 

I  mw 1 *-2*X* Y 

!  w  <  z  >  *  1  -  z 

440 

SUBEND 

line  410 


st.  e epest 


O  n 
3 


450  ' 

460  SUB  B i rec t i on tRewl , I rowl , Cos, S i n >  !  direction  of 

470  T*SQR<Rewl*Rewl  +  Iiriwl*Imwl  >  !  steepest  descent 

480  IF  T>0  THEN  510 
490  Cos*Sm*0 


500  GOTO  530 
510  Cos*-Rewl/T 
520  Sin*ImwlxT 
530  SUBEND 
540  ! 


descent 


A-3 
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55©  SUE  Mul<Xi,Yl,X2,Y2,A,B>  1  21*22 

5*0  R=X1*X2-V1*Y2 

570  £=X 1  * Y2+X2* Y 1 

580  SUBEND 

5S0  » 

800  SUB  Di'AXl,  Y1,X2,Y2,A,B>  !  21/22 

810  T=X2*X2+Y2*Y2 

820  A  =  < X  1  *X2+Y  1  *Y2  >  si 

8  3  0  B  - Y 1  *  X  2  -  X 1  *  Y  2  >  -  T 

840  SUBEND 

650  1 

880  SUB  E  ».  p  X  Y ,  A ,  B>  !  E  X  P  (.  2 ) 

670  T =EXP ( X > 

680  A=T*COS 1  Y> 

690  B=T*SIN<Y) 

700  SUBEND 
710  ! 

720  SUB  Log<  X, Y, A, B>  »  PRINCIPAL  LOG<  2> 

730  R= . 5*L0G  <  X*X+Y*Y ) 

740  IF  X< >0  THEN  770 

750  B= . 5+P I *3GN( Y  > 

760  GOTO  790 
770  B*ATN<Y/X> 

780  IF  X<0  THEN  B  =  B+P I  * U -2* < Y< 0 )  > 

790  SUBEND 
800  ! 

810  SUE  Sqr'X, Y, A,B>  !  PRINCIPAL  SGR<2> 

320  IF  XO0  THEN  860 
830  A=B=SQR< .5*ABS<Y>> 

340  IF  Y< 0  THEN  B— B 

850  GOTO  970 

860  F=SGR<:5QR<X*X+Y*Y)  > 

870  T=. 5*ATN<  Y--X> 

88G  A=F*COS<T> 

890  B=F*S IN<  T  > 

900  IF  X >0  THEN  970 
910  T=R 
920  A=-B 
930  B=T 

940  IF  Y>=0  THEN  970 
950  A=-A 
960  B=-B 
970  SUBEND 
980  ! 

990  SUB  RrgCX, Y, A>  !  PRINCIPAL  ARG'Z> 

1000  IF  x=0  THEN  A=.5*PI*SGN< Y> 

1010  IF  XO0  THEN  fl*RTN<Y/X) 

1O20  IF  X<0  THEN  A*R+PI* ( 1 -2* <  Y<0>  > 

1030  SUBEND 
1040  ! 

1050  SUB  Power < X , Y , R , R , B >  !  PRINCIPAL  POWER  2  R 
1060  F  =  E XP< .  5*R*L0G-:X*X+Y*Y)) 

1070  CALL  RrgCX, Y, T > 

1080  R*F*COS<R*T> 

1090  B=F*SIN<R*T> 

1100  SUBEND 
1118  ! 
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1120  SUB  Siru'X,Y,fi,B>  !  SIN<2> 

lisa  e=exp-y:> 

1140  fl=.5*SlNkX)»(E  +  l-'E) 

1150  IF  ABStYX.l  THEN  1180 
1160  S= .  5*  *:  E- 1  •  E  ) 

1170  GOTO  1200 
1130  S - Y  *  Y 

1190  S= Y* <  1 20*$*  <  20+S  >  )/'120 
1200  B»COS<X>*S 
1210  SUBEND 
1220  ! 

1230  SUB  Cos<X, Y, A, B>  !  COS<Z> 

1240  e=e;:p<y> 

1250  A=.5*C0S<X>*<E+1''E> 

1260  IF  ABSOYX.l  THEN  1290 
1270  S * .  5 * < E - 1  E ) 

1280  GOTO  1310 
1290  S  =  Y  *  Y 

1  300  S«Y*<  1 28  +  S*  <  20+S  >  ')  /  120 

1310  B*-SIN<X>*S 
1320  SUBEND 
1330  ! 

1340  SUB  SinhCX, Y,fl,B>  !  SINH<Z> 

1350  E=EXP<X) 

1360  B=.5*SIN(Y)*<E+1/E) 

1370  IF  ABS <  X  >  <  • 1  THEN  1400 
1380  S=.5*<E-1/E> 

1390  GOTO  1420 
1400  S=X*X 

1410  S=X*  U  20+S*  <  20+S  '»>/ 120 
1420  A=COS  <  Y  )  *S 
1430  SUBEND 
1440  ! 

1450  SUB  CosMX,  Y,  A,  B>  !  COSH<Z> 

1460  E=EXP  <  X  > 

1470  A= . 5*C0S< Y ) *<E+ 1 /E  > 

1480  IF  A B S < X > < . 1  THEN  1510 
1490  S*.5*<E-1/E) 

1500  GOTO  1530 
1510  S=X*X 

1520  S*X*< 120  +  S*<20+S>  >' 120 

1530  B  =  S I N  ( Y  ■•  *  S 
1540  SUBEND 
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Appendix  B 

Steepest  Descent  for  General  Analytic  Function 


Here  we  do  not  force  the  integrand  to  be  of  the  form  in  (I),  but  consider  the 
general  integral  /(  dz  f(z).  Let  f(z)  be  analytic  in  a  region  in  the  complex  plane.  The 
magnitude-squared  value  is 


M  =  |f(z) |2  2  f'  ♦  f2 
r  i 


(B-l) 


The  direction  of  steepest  descent  for  M  is  opposite  to  the  gradient  of  M,  which  is 

VM  =  —  a  +  — -  a  .  (B-2) 

3x  x  3y  y 

But,  from  (B-l), 


3M  -  2(f 

3f 

r 

3x  \r 

3x 

2H.  2(f 

3f 

r 

3y  l  r 

3y 

df . 

+  fi  IT 


3f . 


+  f. 


_i \ 

3y  j 


Now  if  function  f  is  analytic  at  z,  then 

3f 


„  3f.  3f. 

. . ,  ,  r  .i  i 

f  (-Z')  3x  +  X3x  =  3y 


3f 

i 

3y 


(B-3) 


(B-4) 


So  we  can  express 


=  2  Re{f*(z)  f'(z)} 

=  -2  Im{f*(z)  f'(z)} 
dy 


(B  -5) 


Therefore,  the  steepest  descent  direction  for  M  has  components  that  are  the 
negatives  of  (B-5)  or,  equivalently,  are  proportional  to 


-Re 


f(z) 


r  ;  f  ^ 


(B-6) 


The  basic  calculation  for  determination  of  steepest  descent  directions  is  thus  seen  to 
be  , 

fWffz)  =  ~j~  In  f(z). 

Whenf'(zs)  =  0,  we  have  a  saddle  point  of  f  at  z  =  z5.  Near  the  saddle  point, 

f(z)  =  f(z  )  +  if"(z  )(z  -  z  )2  e  £  +  if  .A2  for  small  A  .  (B-7) 

s  2  s  s  o  2  2 


B-l 
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Now  let 


1,  . 2\ /  *  1 


f.  1 *  A*  ) 


K  v2} 


for  small  A 


Then  (B-8)  yields 


C  c  la  ,  10 

f  f  =  ac  ,  A  =  re 
o  2  5 


2  ? 

M  =  |f  |“  +  ar  . cos(a  +  29)  , 


(B-10) 


which  has  two  peaks  and  two  valleys  versus  0  in  a  2n  interval  (for  a  #  0);  this  is 
characteristic  of  a  saddle  point.  The  directions  of  steepest  descent,  0d,  at  z  =  zs 
occur  when 


<  +  20 .  =  *  or  3a  ,  G 


7T  -  Qt  71  -  a 

— 75 —  or  — ^ —  +  11 


(B- 11) 


Notice,  from  (B-9),  that 


'  =  }  =  arg{f  (:s)  =  ars{f"^s^/f (zs)}  •  (B-l  2) 

Nov.  let  us  investigate  the  behavior  of  the  complex  function  f(z)  near  zs,  along  the 
steepest  descent  contours.  For  0d  as  given  by  (B-l  1),  A  in  (B-9)  becomes 


A  =  +.r  exp 


.2  2  -ia 

A  =  -r  e 


and  there  follows,  from  (B-7)  and  (B-9),  for  small  r, 

f(z)  =  f  -  -U2f9e"la  =  f  (l 
o  2  2  o 


on  steepest  descent  contour  near  z( 
That  is,  since  the  term  in  parentheses  is  real  and  positive. 


(B-13) 


(B- 14) 


arg  f(z)  =  arg  f(zs)  on  steepest  descent  contour  near  zs*  (B-15) 
More  generally,  it  can  be  shown  that 
arg  f(z)  =  arg  f(zs) 


everywhere  on  steepest  descent  contour  through  zg. 


(B-l  6) 


Now  let  us  apply  these  general  results  to  the  special  case  where  f(z)  =  exp(w(z)). 
Then  f '(z)  =  w  '(z)  exp(w(z)),  and  f '(z)  =  0  when  w  '(z)  =  0.  Thus  saddle  point 
locations  are  as  usually  stated.  Also,  as  needed  in  (B-6), 


B-2 
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f '  (7)  1 

frir  =  *  (z)  *  (B-17) 

and  then  (B-6)  agrees  with  (9).  Furthermore,  since 

f"(z)  =  [w"(z)  +  w,2(z)]  exp(w(z) )  , 

f,,(zs)  =  w"( zsD  exp(w(zs)) 

then  (B-12)  yields 

a  =  arg  w  (zs) 

When  (B-19)  is  used  in  (B- 11),  the  steepest  descent  directions  corroborate  reference 
1 ,  (7. 1 .8)  and  (7. 1 . 19).  Finally,  (B-16)  yields 

arg  f(z)  =  arg{ exp (w(z))}  = arg{exp(u  +  iv)}  =  v  ,  (B-20) 

meaning  that  v  is  constant  on  steepest  descent  contours;  this  agrees  with  reference  1 , 
Lemma  7.1. 

What  this  all  demonstrates  is  that,  for  a  general  given  integrand  f(z),  we  can  let 
w(z)  =  I  n  f(z)  and  apply  our  usual  techniques  on  w(z).  This  procedure  was  adopted 
in  the  Gaussian  exponent  example  in  the  main  text  of  this  report. 


(B-18) 


(B-19) 
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